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Abstract 

The gauge usually adopted for extracting the reduced Hamiltonian of a thin 
spherical shell of matter in general relativity, becomes singular when dealing with 
two or more intersecting shells. We introduce here a more general class of gauges 
which is apt for dealing with intersecting shells. As an application we give the 
hamiltonian treatment of two intersecting shells, both massive and massless. Such 
a formulation is applied to the computation of the semiclassical tunneling probability 
of two shells. The probability for the emission of two shells is simply the product of 
the separate probabilities thus showing no correlation in the emission probabilities 
in this model. 
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1 Introduction 



A lot of work has been done in the subject of thin spherical shells of matter in general 
relativity [U E] and several applications given [31 HI [5]. An interesting application of the 
mechanics of thin shells has been the semiclassical treatment of the black hole radiation 
[HIEIIH]- An important role in such a treatment is the choice of the gauge. In the original 
treatment [6] a gauge was adopted in which the radial component of the metric is equal to 
the radial coordinate except for an arbitrarily small region around the shell. Such a gauge 
choice was put in mathematically clear terms in [9l [10] where the reduced canonical 
momentum is extracted through a well defined limiting process. However such a limit 
gauge becomes singular when two or more shells are present and intersect. In addition 
also in the simple instance in which only one shell is present, it would be nice to have a 
procedure in which no limit process is present as the gauge choice should be completely 
arbitrary provided it is free of coordinate singularities. Moreover even in the one shell 
case the procedure for extracting the canonical momentum is usually considered as a very 
complicated one. Here we shall give a treatment which greatly simplifies the derivation 
of the reduced action, does not require any limiting procedure and can be applied for the 
treatment of two or more shells which intersect. 

All the problem is to derive the reduced action, i.e. an action in terms of the shell 
coordinates and some appropriate conjugate momenta. We shall keep the formalism for 
the two shell close as possible to the one shell treatment. 

We shall perform all the treatment for a massive dust shell and the simpler case of a null 
shell can be derived as a particular case. As an application of the developed formalism 
we rederive the well known Dray-'t Hooft and Redmount relations for the intersection of 
light-like shells PHQ2]. 

The main motivation which originated the works 0, [8j [13] is the computation of the 
semiclassical tunneling amplitude, which is related to the black hole radiation in which 
one takes into account the effect of energy loss by the black hole. It is remarkable that 
such a simple model reproduces all the correct features of the Hawking radiation, giving 
also some corrections due to energy conservation. For approaches more kinematical in 
nature see e.g. [H]. The adoption of more general gauges allows the dynamical treatment 
of two intersecting shells without encountering singularities. The formalism developed 
here is applied to the computation of the semiclassical emission probability of two shells; 
it was in fact suggested [6] that correlations could show up in the multiple shell emission. 
We find however that such a model gives no correlation among the probabilities of the 
two emitted shells. 
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The paper is structured as follows. 

In Section 2 we lay down the formalism by exploiting some peculiar properties of a function 
F strictly related to the momenta canonically conjugate to the metric functions. It is then 
just a simple matter of partial integrations to extract the reduced canonically conjugate 
momentum without employing any limit process. This is done in Section 3. A new term 
containing the time derivative of the mass of the remnant black hole appears; if such a 
mass is considered as a datum of the problem the result agrees with those of the original 
massless case of Kraus and Wilczek and with the massive case result obtained through 
a limiting process in [9]. Then we discuss the derivation of the equations of motion. In 
[U] and [U] the variational principle was applied by varying the total mass of the system 
which we shall denote by H. In [S] the attitude was adopted of keeping the total mass 
of the system fixed and varying instead the mass of the remnant black hole. Here we 
show that both procedures can be applied to obtain the equations of motion; depending 
however on the choice of the gauge one procedure is far more complicated than the other 
and we give them both. 

In Section 4 we discuss more general gauge choices and derive the equation of motion in 
the inner gauge. We shall not consider in the present paper complex gauges or complex 
gauge transformations [15J. 

In Section 5 we discuss the analytic properties of the conjugate momentum; this is of 
interest because the imaginary part of the conjugate momentum is responsible for the 
tunneling amplitude and in determining the Hawking temperature both via a simple me- 
chanical model or more precisely by working out the semiclassical wave functions on which 
to expand the matter quantum field (TU]. We remark that the result for the tunneling 
probability is independent of the mass of the shell but depends only on the initial and 
final energy of the black hole. 

In Section 6 we extend the treatment to two shells in the outer gauge writing down the 
explicit expression of the two reduced canonical momenta. 

In Section 7 we derive the equations of motion for both shells from the reduced two shell 
action. For simplicity this is done for the massless case; the general treatment for massive 
shells is given in Appendix C. 

In Section 8 the developed formalism is applied to give a very simple treatment of two 
shells which intersect. In the massless case we rederive the well known relations of Dray 
and 't Hooft [H| and Redmount p2]. 

In Section 9 we consider the problem of computing the tunneling probability for the emis- 
sion of two shells which in the process can intersect. To this end one has to compute the 
imaginary part of the action along the analytically continued solution of the equations of 
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motion. In this connection a helpful integrability result is proven which allows to compute 
the action along a specially chosen trajectory on the reduced coordinate space. Such a 
result allows on one hand to prove the independence of the result from the deformation 
defining the gauge and on the other hand it allows to compute explicitly such imaginary 
part. The final outcome is that in all instances (massive or massless shells) the result 
again depends only on the initial and final values of the masses of the black hole and the 
expression coincides with the one obtained in the one shell case. The interest in studying 
the two shell system was pointed out in [6] in order to investigate possible correlations 
among the emitted shells. Here we simply find that the two shells, even if they interact 
with an exchange of masses, are emitted with a probability which is simply the product 
of the probabilities for single shell emissions and thus that the model does not predict 
any correlation among the emitted shells. 
In Section 10 we summarize the obtained results. 



2 The action 

As usual we write the metric for a spherically symmetric configuration in the ADM form 

mm 

ds 2 = -N 2 dt 2 + L 2 (dr + N r dt) 2 + R 2 dU 2 . (1) 

We shall work on a finite region of space time (ti,tf) x ( r o 5 r m) • On the two initial and 
final surfaces we give the intrinsic metric by specifying R(r,ti) and L(r,ti) and similarly 
R(r,tf) and L(r,tf). 

The complete action in hamiltonian form, boundary terms included is [T7| [3], [9] 

rtj j>r m _ _ pt f NRR' T ' m 

S = S shell + dt dr(n L L + ir R R - NH t - N r H r ) + / dt (-N r n L L + — — ) 

Jti Jr Jti L r 

(2) 

where 

/"*/ 

S s heli = dtpr. (3) 
Jti 

Tit and TC r are the constraints which are reported in Appendix A, f is the shell position 
and p its canonical conjugate momentum. Action (j2J) is immediately generalized to a finite 
number of shells. S s hdi as given by eq.flB]) refers to a dust shell even though generalizations 
to more complicated equations of state have been considered [21 [31 [18]. 
Varying the action w.r.t. iV and iV r gives the vanishing of the constraints while the 
variations w.r.t. R, L, 7r#, ttl give the equations of motion of the gravitational field [H 
|9] which for completeness are reported in Appendix A. The functions R, L, N, N r are 
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continuous in r while R', L', N', N r ' , tt^, itr can have finite discontinuities at the shell 
position [5]. In [U] it was proven that the equation of motion of a massive shell cannot 
be obtained from a true variational procedure, in the sense that the obtained expression 
for p is discontinuous at r = f. In the same paper it was remarked that for consistency 
the average value of the r.h.s. of the equation for p has to be taken. For the reader's 
convenience we give in Appendix A the explicit proof that the equations of motion of 
matter i.e. f and p can be deduced from the already obtained equations of motion for the 
gravitational field combined with the constraints. The equation for f does not pose any 
problem while for the equation for p one deduces algebraically that the r.h.s. contains 
the average of the derivatives of L, N, N r across the shell. In Appendix A we also show 
directly that such discontinuity is absent in the massless case |lUj . 

Already in P E] it was pointed out that in the region free of matter, as a consequence of 
the two constraints the quantity M. 

. . 7r? R RR' 2 , . 

M = — H r- (4) 

2R 2 2L 2 K J 

is constant in r and this allows to solve for the momenta [21 E] 




k l = //« (-) -1 + — =RW (5) 

L[(R/L)(R'/L)' + {R'/Lf - 1 + M/R] ff ,. 

7[R = w ^ 

In the one shell problem we shall call the value of such A4, M for r < f and H for r > r. 
The function 



\R'\ 2 _ 2M IK I f R'Y . 2M 



hl ^t) - 1 + ^r +RR '^[T-i{T) -i+- s -l-*/'w P) 

has the property of generating the conjugate momenta as follows 

OF 

« ^ 3L (»» 
6£_dF__d_dF_ 
nR ~ 6R~ dR~ IfrdR/' ( ' 

The total derivative 9 ^ = R'f'(R) of the arbitrary function f(R) does not contribute 
to the momenta. The function F will play a major role in the subsequent developments. 
A large freedom is left in the choice of the gauge. With regard to L we shall adopt the 
usual gauge L = 1. It will be useful in the following developments to choose the arbitrary 
function f(R) such that F = for R' — 1. Such a requirement fixes f'(R) uniquely. 

F(R, R', M) = RW(R, R', M) + RR'(C(R, R' , M) - B(R, M)) (10) 
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where 



W(R,R',M) = \/R' 2 -1 + ^-; C(R,R',M) =log(R' - W(R,R',M)) (11) 



and 



B(R,M) = f-^- + log(l - f-^). (12) 
The function F(R,R',Ai) has the following useful properties 

dF(R,R',M) dC _ 1 

— m — - R{c - b), M ,-- w (13) 

Other properties of F will be written when needed. 

In the following section we shall choose R = r for r > f and also R = r for r < f — I so that 
F vanishes identically outside the interval f — I < r < f. We shall call this class of gauges 
"outer gauges". In Sect. H]we shall consider more general gauges e.g. the gauge R = r for 
r < f and also R = r for r > f + I which we shall call "inner gauges" . However contrary 
to what is done in [HI H] we will not take any limit I —>■ and prove that the results are 
independent of the deformation of R in the region f — / < r < f (or f < r < r + I for 
the inner gauges) provided R' satisfies the constraint at r = f . We shall call these regions 
(f — /, f) for the outer gauge and (f, f + I) for the inner gauge, the deformation regions. 
The variation of S which produces the equations of motion has to be taken, as it is well 
known, by keeping the metric and in particular R and L fixed at the boundaries. The 
variation of the boundary terms gives 

- N r (r m )67T L (r m ) + N r (r )57r L (r ). (14) 

The N, N r can be obtained from the two equations of motion for the gravitational field 

= N llk - if] + W; R=-N^ + N r R'. (15) 

Using these it is easily proved that for r outside the deformation region (f — l,r) we have 
iV r = N \ — for r > r, N = const and N r = N\f— for r < f — L N = const where the 
two constants as a rule differ. Thus the variation of the boundary term is 

- N(r m )5H + N(r )5M. (16) 

In the next section we shall connect N{r m ) with N(r ) being N(r) not constant in the 
deformation region. 
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3 The one shell effective action in the outer gauge 



As outlined in the previous section we shall choose 

V(t) 



R(r,t) = r + 



V(t) 



p(r' — r(t))dr' = r H — g(r — fit)) 

r 



m Jo 

having p support (—1,0), p(0) = 1 and smooth in — / and 



(17) 



f 

/ p(r)dr = 0. 



(18) 



As a consequence the deformation g(r) has support in (—1,0) and g'(0 — e) = 1. Such 
R satisfies the discontinuity requirements at r = f which are imposed by the constraints. 
In fact the constraints (see Appendix A) impose the following discontinuity relations at 
f (we recall that we chose L = 1) 



AR' = - 



V 
R ] 



V = a/]3 2 + m 2 



and 



An L = -p. 

In the outer gauge the bulk gravitational action becomes 



I g dt 



(19) 
(20) 

(21) 



where, keeping in mind that L = 1 



/ 



(n L L + n R R)dr 
*® dF 



ro 



irjiRdr 



r dt 

d_ r m 

dt 



dr - M(t) 



' 'OF d OF . 
{ dR~d-rdR )R dr 



ro 



dM 



OR' 



ro 
f(t) 

ro 



Fdr 



ro 



rr{ 

- M(t) / 

Jr 



Ht) Qp 

dM 



dr — 



dF ■' 



(22) 



f(t)-e 



where we used the fact that F vanishes at r = ro. 

Adding l s he\i = W and neglecting the total time derivative which does not contribute to 
the equations of motion, we obtain for the reduced action in the outer gauge 



Pet- M(t) 



ro 



^ dr+ (-N r n L + NRR')\'™ \dt 



(23) 
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where using (j!9|20p 

1 dF 

+ p 



Pc = -F(r(t) - «) " ^ 



f(t)-S 



(24) 



= v^-y^-flog( f + ^ 2 + m2 ~^_- v/ ^ 

y f - V2Mf 

A few comments are in order: 1) No limit / — > is necessary for obtaining p c of eq.(T5] 
which holds for any deformation g. 2) The M(t) term is important, as we shall see, if we 
consider the variational problem in which M is varied. On the other hand if we consider 
M as a datum of the problem and vary H the contribution M(t) is absent. 3) p in eq. (1241) 
is a function of f, H and M as given by the discontinuity relation ( 1201) equivalent to the 
implicit equation 

H-M = V + — -m— (25) 
2r V r 

We discuss now these issues in more detail. Let us consider at first M as a datum of 
the problem and vary H . As shown in Appendix A in order to be consistent with the 
gravitational equations M has to be constant in time. This is the situation examined in 
[9] where the expression (1241 for p c was derived by a limit process in which I — > 0. From 
eq. (1231) we see that the equation of motion for f is given by 

*W ~ N{Tm) = (26) 

where N(r m ) can be replaced by N(f) being N(r) in the outer gauge constant for r > f. 
The computation of eq. (121)]) keeping in mind the implicit definition of p eq. (1231) gives the 
correct equation of motion for the massive shell 

f = LN{f) - N r (r) =UL- N(r). (27) 

The outline of the calculation is done in Appendix B. 

Alternatively one can consider H = const as a datum of the problem and vary M(t). We 
remark that as shown in Appendix A the datum M or H is consistent with the gravita- 
tional equations only if H and M are constant in time. Nevertheless in the variational 
problem H and M have to be considered as functions of time, because the constraints tell 
us only that M and H are constant in r. Only after deriving the equation of motion one 
can insert the consequences of the gravitational equations of motion. 
The variation of M(t) in the outer gauge is a far more complicated procedure, due to the 
presence of M(t), but gives rise to the same result obtained by varying H and keeping M 



7 



fixed. As this will be useful to understand the two shell reduced dynamics to be developed 
in Sect JH] we go into it with some detail. In this case the M term plays a major role; in 
fact the equation of motion now takes the form 



; dp c d 
~dM + di 



dr + N (r )777T( r o) 



''0 



dM 



dM 







(2* 



where due to the vanishing of g(x) for x < — I we have vr^(r ) = i/2Mr . N r (r) on the 
other hand is obtained by solving the two coupled equations (JT5|) with the condition that 
for r > f , N T (r) equals having normalized N = 1 for r > f. One easily finds that 

for r < f 



N r (r) 

Taking into account that 



W 



d 



R^dr + 
dM 



d 2 F . d 2 F 

+ 



V2H~7 



V2Hf + p 



dr y dR'dM J dMdR dM 



we have 



and eq. 



d r m dF 
df 

becomes 



dM 



dr 



r R^dr + 



dM 



: dF 



d 2 F 



dM dMdR 



R 



a dp c 



dM 



+ 



x dF ■ d 2 F 
T dM + dR'dM 



+ 



V2W7 



0. 



From the expression for p c given by the first line of eq. 

dp V dR 



dM W dM 



V2Hi 



\ / 2Hr + p 

we obtain 

- = 



(29) 

(30) 
(31) 

(32) 



V2H: 



P 



where W(r — e) 
Using 



'R! 2 {f + 2M/R(r) 



dM 



p. 
v 



(33) 

+ I and R'(f - e) = 1 + V/f. 

(34) 

is necessary 



we obtain eq.l jzTj) again. We remark once more that no limit process / 
for all these developments. 

To summarize, in the present section we have derived the reduced action for the one shell 
problem in the outer gauge with an arbitrary deformation. One can vary H (the exterior 
ADM mass) considering the interior mass as given, or one can vary the interior mass M 
considering the exterior mass H as given, or even one can vary both M and H always 
obtaining the correct equations of motion. Whenever M is varied the M term in eq. (1281) 
plays a crucial role. All the results do not depend on the deformation g. 
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4 More general gauges 



It is of interest to examine more general gauges given by eq. (117ft where g does not neces- 
sarily vanish for positive argument, i.e. we can consider g(x) with g(x) = for \x\ > I, 
g(0) = and g'(+0) —g'(—0) = —1, thus again satisfying the constraint (fl9"j) . In this case 
F(R, R f , H) does not vanish identically for r > f and the bulk action (j2ip is given by the 
time integral of 



OF ■ 



r(t)+e 

(35) 

f(t)s 



d F m ■ r (t) dF ■ r Vm OF 

The p c is easily computed using eq. (1351) and one gets immediately the general form for 
the canonical momentum p c 

p c = f(A£ - AB) (36) 
where AC = £(f + e) — £(f — e) and similarly for AB and the reduced action becomes 



f f . . r f W dF ■ f Tm dF 



S ----- / ( p c r - M(t) J —dr - Hit) —dr + (-N r n L + NRR')\ ™ \ dt. 

(37) 

We shall call inner gauge the one characterized by g(x) =0 for x < 0. 

Due to the similarity with the treatment of Sectj3] we shall go through rather quickly. 

Now F vanishes identically for r < f and the action takes the form 

f tf ( ■ ■ ■ F m dF \ 
S = jf dt (p»f -H J —dr + (-N r 7T L + NRR')\ r ™j (38) 



where p % c is given by 

i [„ R OF 

v — F + 

Pc [ f OR' 

whose explicit value is 



+ p (39) 

f+e 



pi = V2Mr - V2Hf - f tog ( — ^ 2H r ] (40) 

\f -V+P-V2M? J 

and p, again determined by the discontinuity equation ( 120|) . is given by the implicit 
equation 

H-M = V-£-pM (4!) 



which is different from eq. (1251) . Now the simple procedure is the one in which one varies 
M keeping if as a fixed datum of the problem. This time the solution of the system 
eq. (1151) gives simply iV = const and N r = Na/— for r < r < f. 



9 



5 The analytic properties of p c 



We saw that in the outer gauge 



p c = V2M r — \J2H f — f log 



f + V-p- \/2Hr 



(42) 



f - \/2Mr 



The solution of eq. (1251) for p is 



P 



r 




(43) 



where 



.4 



H - M m 2 



(44) 



2f 2 ' 



If we want p to describe an outgoing shell we must choose the plus sign in front of the 
square root. Moreover the shell reaches r = +oo only if H — M > m as expected. 
The logarithm has branch points at zero and infinity and thus we must investigate for 
which values or f such values are reached. At f = 2H, p has a simple pole with positive 
residue; then the numerator goes to zero and below 2H it becomes 



where here V is the absolute value of the square root. Expression (|45p is negative irrespec- 
tive of the sign of p and stays so because p is no longer singular. In order to compute the 
tunneling amplitude, below f = 2H we have to use the prescription |8J f—2H f—2H—is 
and as a consequence the p c below 2H acquires the imaginary part inr. Below f = 2M 
the denominator of the argument of the logarithm in eq. (1421) becomes negative so that 
the argument of the logarithm reverts to positive values. Thus the "classically forbidden" 
region is 2M < f < 2H independent of m and of the deformation g and the integral of 
the imaginary part of p c for any deformation g is 



with uj = H — M which is the Parikh-Wilczek result [Hj. A more profound way to relate 
the result for p c to the formula for the Hawking radiation is to use (]24l to compute 
semiclassically the modes on which to expand the quantum field, and then proceed as 
usual by means of the Bogoliubov transformation. This was done in j6l [16] ■ A more 
direct particle like interpretation of (I46p was given e.g. in |19j . 




(45) 




(46) 
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Similarly one can discuss the analytic properties of the conjugate momentum p l c in the 
inner gauge. We have 



1 — — 



p- = V2Mf - V2H? - f log ¥-= . (47) 

I _ V _ 2M _|_ " 



This time the solution of eq. (|4ip gives for p 



p J™A±JA 2 -(1- 2M) ^ 



2M 



f 1 - 

r 



(48) 



where now 

t H — M m 2 

A =~i- + 2Fr (49) 
Notice that for m = we have p c = p l c . To describe an outgoing shell the square root in 
(1481) has to be taken with the positive sign. 
All the point is the discussion of the sign of the term 



V 2M p „ . / Q/ 2H 

+ £ = + e)- J Ri2(f + £ )-i + — (50) 



V T* IT V 

where 



2iW T) In/ \ 

'i? /2 (r + e) - 1 + — . (51) 



For m — at r — 2H eq. (|51j) is negative so that at f = 2if 



^ /2M p . / „. 2H 

— + t = i?' (r + e) + \/i?' 2 r + e - 1 + — 52 

is positive, being in eq. (!5*2l the square root on the r.h.s. understood as the positive 
determination of the square root. The same happens for m ^ provided m < H — M 
which is the condition for the shell to be able to reach f = +oo. For f < 2H as —1 + ^ > 

r 

such a term stays positive irrespective of the sign of R', until R' diverges. This happens 
at 2M where p given by eq. fj4*H|) has a simple pole with positive residue. Thus below 2M 
eq. (!5*2l) reverts to 



V 2M p ,, . / n/ 2H 
1 - - - \ — + ~ = i?'(r + £ - \/ J R /2 (r + e - 1 + — 53 

which is negative irrespective of the sign of R'(f + e). The conclusion is that p l c takes the 
imaginary part iirr in the interval 2M, 2H as it happens for p c . 
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6 The two shell reduced action 



From now on we shall work in the outer gauge. We denote with fi and f 2 the coordinates 
of the first and second shell fx < f 2 - The value of M. for r < fx will be denoted by M, 
for fx < r < f 2 by M and for r > f 2 will be denoted by H as before. In extending 
the treatment to two interacting shells we shall keep the formalism as close as possible 
to the one developed in Sect. EJ The most relevant difference is that in any gauge the 
intermediate mass Md and the total mass H always intervene dynamically. We shall 
consider the mass M datum of the problem. 
For the metric component R we shall use 

V 2 Vx 

R(r) = r + v 2 g(r - f 2 ) + vxh(r - fx); v 2 = v 1 = (54) 

R{r 2 ) R(rx) 

where h(x) has the same properties of g(x) described in Sect. |3] (actually we could use 
the same function). Both g and h vanish for positive argument and thus we are working 
in an outer gauge according to the definition of Sect. [3] 
The action is given by the time integral of 

I = Vi^i + ftih + J dm R R + b.t. (55) 

Breaking the integration range from to fx and from fx to f 2 and using eq.Q for 7Tr 
and the same technique as used for the one shell case, we reach the expression 

P,h + Pxh + Jf J o Fdr - M I —dr - f 2 F(f 2 - e) - (R— j (f 2 - e) + 

(dF \ 
Rq^) (fx) (56) 

d f^ 2 H 2 dF ( dF \ 

= M + * L Fir ~ ^ L am dr + p ° 2 h + h AF(h) + A { ti aW) ih) 

where p° c2 ls given by eq.( l24"l) with M replaced by Mo, and A stays for the jump across 
the discontinuity at fx- We recall that F depends only on R,R' and A4 and the partial 
derivative w.r.t. these variables will have the usual meaning. For the other quantities 
appearing in the calculations we recall that the independent variables are M, M , H, fx, 
f 2 and when taking partial derivatives we shall consider the remaining variables as fixed. 
In eq.( |56l) we have, using eq.( |T3i) . denoting with A the discontinuity across fx and with 
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the bar the average at ri 



dF 

hpi + fi AF + A ( — R ) = (57) 



where 



dF ■ — dF 
= r^ l+ r 1 AF + —AR + RA— = 

h[R'R{AC - AB) + AR'R(£ - B)\ + AR R{£ - B) + RR(A£ - AB) = 

dT ■ dT 
= hPd + hPc2 + H(R(h) - ri)gjjV + M (R(h) - ^)g]^ V 

T = logv 2 ; V = R(A£ - AB); p cl = R\h + e)V; (58) 

dT d 

Pc2 = -{R'{h + e) - l)V + (R(h) - h)-&rT> = -jr{R{h) ~ h)V (59) 
having used 

AR = -hAR' (60) 

and 

R = -hj - h(R' - l) + {R{h) - ri)— . (61) 

Summing up the reduced action for the two shell system is given, boundary terms included, 
by the time integral of 

dT ■ dT 

d r 2 ■ r 2 dF 
+ Jt I Fdr ~ M " L m dr + ( ~ N '* L + NRRX: (62) 

where p c2 = p° c2 + p C 2 with jP c2 given by eq. fT2^I) with M replaced by M and p c2 by eq. 
( |59|) . With regard to eq. fl6"21) we notice that irrespective of the gauge used both terms in 
Mo and H appear in the action. Moreover p c i and p c2 depend both on f\ and f 2 . 



7 The two shell equations of motion 

From the reduced action f)62p we can derive the equations of motion for the two shells. 
This is of some importance in order to show the consistency of the scheme. We recall 
that action (1621) has been derived in the outer gauge i.e. R(r) = r for r > f 2 (f 2 > f\). 
While in the one shell problem formulated in the outer gauge H does not appear, in the 
two shell problem it is always present. Again we consider M = const, as a datum of the 
problem. 



13 



In the variational procedure we can vary f 1; f 2 , H and M independently. We shall start 
varying H but keeping all other parameters fixed. For the sake of simplicity we shall deal 
here with the massless case m% = m 2 = 0. In Appendix C we give the general derivation 
for mi and m 2 different from zero. The most important fact that occurs when we vary 
H keeping Mo fixed is that the terms proportional to f\ cancel in the variation. The 
simplifying feature of the massless case is that AC = so that in eq. (!581) V = —RAB, 
being B function only of R and M. and not of R' and thus only of R and M being M 
a datum of the problem. The coefficient of r\, taking into account the following relation, 
easily derived from ( fl3l) . 

dF + R'R dE - 1 ffi* 
dM +RR dM~W^R (63) 

is found proportional to 

^(R'(h + e)V) - {R!{h + £ )~ 1 )§j V - WO - VjfiM^ 1 + £) (64) 
where we used 

^pl = f lR ' {h +e)- h{R\h + e) - 1) + ^Wi) - h) (65) 

and we took into account that T does not depend on fi and that on the equations of 
motion H = M = 0. Now employing the relations 

&R(ri) &T _ dR'ih + e) dT .... . 

ir=s«"»-'* — r = m Wri+ ^-^ (66) 

we see that expression (|64|) vanishes. Thus we are left only with the f 2 terms. We know 
already that the boundary term and the p° c2 term give the correct equation of motion and 
thus we have simply to prove that the f 2 terms originating from 

-j t ^-^m v] (67) 

, dp c2 . 
cancel — — i.e. 
on 

-^[-{R'ih + e)- 1)V + (R(h) - h)§-V] - ^P(ri) - h)j^V] = 0. (68) 

Using relations ( !66l) we have that eq.( l68l) is satisfied. Such a result is expected as the 
exterior shell parameterized by f 2 moves irrespective of the dynamics which develops at 
lower values of r until fi crosses f 2 . 

Now we vary M keeping H fixed. We have no boundary term contribution because also 

M is constant and we find the equation 

■ d Pcl , dp c2 d . dT . ■ dF(f 2 -e) ■ dF(h+e) 
ri dM + r2 dM ~ Jt [{R{ri) ~ ri) dM Vl + r2 ^M ri dM + 



1 R dM dr + R dM^dR 



ri 



= 0. (69) 

fi+e 
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Using expression (|29|) for N r (r), relation (T151) for N(r) and relation (163|) we obtain with 



[R'ih +£)- W(h + e)]-\h + AT r (fi) - N(h)) 



'«9M n 



'• OF d 2 F ■ 
r2 aMo + dM dR> 



+ a 



-iR 



-[R'ih + e) - W(r. ■ :)] ' — (r • :] 



+r 2 



d 2 F 



dM dM dR 



R 



ri+e 



(70) 



The first line is the equation of motion for fx, the second line vanishes being the equation 
of motion for f 2 (see eq. (1321) of the one shell problem) while exploiting the relation 



d 2 F n 8B 

+ R- 



OMdR' 



(71) 



dM W(W - R>) 

also derived from f|T3l and substituting into f 2 the value given by the equations of motion 
derived previously we find that the third line simplifies with the fourth. Thus we are left 
with the relation 

fi + iV r (fi) — N(ri) — (72) 
which is the correct equation of motion for the interior shell. 



8 Exchange relations 

In the present section we consider in the above developed formalism the intersection of 
two shells during their motion. In the massless case we shall rederive the well known 
relations of Dray and 't Hooft [TT] and Redmount [12] between the masses characteristic 
of the three regions before and after the collision. We shall denote by to the instant of 
collision, fo the coordinate of the collision. It is well known that some hypothesis has 
to be done on the dynamics of the collision which should tell us which are the masses of 
the two shells after the collision. Once these are given the problem is reduced to that of 
a two particle collision in special relativity. The formulas which we develop below give 
the intermediate mass M' Q after the collision. The simplest assumption is that of the 
"transparent crossing" i.e. after the collision shell 1 carries along with unchanged mass 
mi and the same for shell 2. A further simplification occurs in the massless case i.e. 
when the two massless shells go over to two massless shells. We develop first the general 
formalism which applies also when we have a change in the masses during the collision 
and then specialize to particular situations. 
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Just before the collision i.e. at t = t — e assuming that fi(t — e) < r 2 (t — e) we have 
for the momentum 7T^ 



n L (f 1 -0,t o -s) = n L (f 1 + 0, t - e) + p x = (73) 
= 7r L (f 2 - 0, t - e) + pi = 7r L (r 2 + 0, t - e) + Pi + Vi 

and after the collision i.e. at t = t + e we have with r 2 (t + e) < fi(t + e), 

7TL(r 2 -0,t + £) =rc L {f 2 + 0,t + e)+p' 2 = (74) 
TrzXn -0,t + e) + p 2 = 7r L (fi + 0, t + e) + Pi + p' 2 

and at the collision f 2 = fi = fp. The main point in treating the crossing is to realize 
that the sum V\ + V 2 has to be continuous in time as it represents the discontinuity of R' 
at the time of crossing. In fact just before the crossing we have for r < f , choosing for 
simplicity in (!54l) g(x) = h(x), 

Mr) =r+ ^ Vl t V ^ g(r - r ) (75) 
while immediately after we have 

(VI + V') 

R ( r ) = r+ K *Z 2) g(r-r ). (76) 
r o 

If Vi + V 2 7^ V{ + V 2 we would have a discontinuous evolution of the metric for r < r which 
is incompatible with the equations of motion of the gravitational field. As M is unchanged 
during the time evolution Vi + V 2 = V{ + V 2 implies that 7Tx,(ri — 0, to— e) = vr^(f 2 — 0, t +e) 
which combined with the fact that t^l{j'2 + 0, to — e) = ttl(^i + 0, to + e) = \/2Hr Q gives 
the further relation pi + p 2 = p[ + p' 2 . Thus we have the same kinematics as in the two 
particle collision in special relativity. In the case of conservation of the masses mi and 
m 2 we have the same kinematics as that of an elastic collision in 1 + 1 dimensions. A 
discussion of special cases with massive shells using a different formalism has been given 
in [20]. 

The relation between n^fi — 0, t — s) and n L (f 2 + 0, t — s) 



^i + V 2 \ 2 1 2M 




1 + — = V2Hf +p 1 +p 2 (77) 
?"o J ro 

gives the equation 

m\ + ml + 2(V X V 2 - pip 2 ) + 2(Vl + V 2 )r + 2Mf = 2Hf + 2(p x + p 2 ) y / 2Hr~ (78) 
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where M is given by 



H-M^K + ^-n^. (79) 



From the knowledge of p\ one derives M' Q from 



For "transparent" crossing i.e. m'- = rrij and p'- = pj, M' Q is given by (IHUj) with = p\. 
The massless case is most easily treated. In this case in order for the shells to intersect 
they must move in opposite directions, the exterior one with negative and the interior 
one with positive velocity. The conservation of V\ + V2 and p\ + p2 in the massless case 
has two solutions which are physically equivalent i.e. p\ = pi, p' 2 = p 2 - The two shells 
intersect only if p 2 < and f>\ > so that V2 = —p2 and V\ — Pi- From eq. (J79l) we have 

H-M 

P2 = 1= 81 

and from p[ = p\ 

H — ML 

Pi = /==■ 82 

1 — — 

Eq. flTS^) now becomes 

Hf + V2~H^o(pi + P2) = f {pi - P2) + Mr - 2pip 2 (83) 
and substituting here eqs. fl81|82p we obtain 

Hr + Mf - 2HM = M r - 2M M' + M^f (84) 
which is the well known Dray- 't Hooft and Redmount relation [TTl [T2] . 

9 Integrability of the form p c idfi + £> C 2<^2 an d the two 
shell tunneling probability 

We are interested in computing the action for the two shell system i.e. the integral 

/ dt(p cl fi + p c2 f c2 ) = / (pci dfx + p c2 rff 2 ) (85) 

•Jti Jru,r2i 

on the solution of the equations of motion. This is of interest in the computation of the 
semiclassical wave function and the tunneling probability in the two shell case. We shall 
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prove explicitly that the form p c \dr\ + p C 2dr 2 is closed i.e. integrable. Such a result will 
be very useful in the actual computation and in showing the independence of the results 
of the deformation g. In fact we can reduce the computation to the integral on a simple 
path on which the two momenta p c \ and p c2 take a simpler form. 

The integrability result is similar to a theorem of analytical mechanics [211 122] stating 
that for a system with two degrees of freedom, in presence of a constant of motion the 
form pidqi + P2dq 2 where the pj are expressed in terms of the energy and of the constant 
of motion, is a closed form. Here the setting is somewhat different as we are dealing with 
an effective action with two degrees of freedom, arising from the action of a constrained 
system. We recall that for f\ < f 2 M Q like H is a constant of motion in virtue of the 
equations of motion of the gravitational field. The p C j are functions of H, M ,f 1 ,f 2 in 
addition to the fixed datum of the problem M. The effective action takes the form 

dt(p c A + pJ-2 + PhH + p Mo M ) + b.t. (86) 

where b.t. is the boundary term (see eq.([2D) which depends only on H,N r and not on 
fj. The value of N r is supplied by the solution of the gravitational equations of motion. 
Thus varying w.r.t. f\ we have 

n dp cl ■ dp cl x dp c2 dp H tt dp m ,v {Qr? , 

= ^ +r x —- +r 2 —- + —-H + ——M . (87) 

at ori ori ori or 1 

We show in Appendix A that the constraints combined with the equations for the gravi- 
tational field impose = H = M and thus we have 

_ = 0. (88) 
<9f 2 df\ 

The meaning of the procedure is that the consistency of the variational principle imposes 
eq. (1881) . On the other hand eq.( l88i) can be verified also from the explicit expression of p c \ 
and p C 2 given in Sect. [6j If the crossing of the shells occurs outside the region 2M, 2H we 
can choose the path as to keep the two deformations non overlapping in such a region; 
then p c i takes the simple form (|2~^|) with H substituted by M and p c2 again the form 
with M substituted by M and thus one obtains for the imaginary part of integral ( 1851) 



i ((2M ) 2 - (2M) 2 + {2Hf - (2M ) 2 ) = i {{2Hf - (2M) 2 ) (89) 

i.e. the two shells are emitted independently. If the crossing occurs at f with 2M < 
f < 2H, with e.g. r\ < r 2 before the crossing, we choose the path f 1 = f 2 — e before 
the crossing and r 2 = fi — e after the crossing. For clearness sake we examine at first the 
problem for the crossing of a null shell with a massive shell, even when the mass of the 
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massive shell changes. In order to reduce the integration path to the described one, one 
must first, given the two initial points with R(fu) < 2M and R(f 2 i) < 2M, bring them 
together ( fi will denote the position of the massless shell). For fu < f 2i this is done by 
integrating along the line f 2 = const with f\ varying from fu to f 2 % — e. The contribution 
is 



Pddh 



R'(f x + e) V dh- 



(90) 



'<'i! 



But R'ifi) is real and 



V 



-R{h) 




2M 

rW) 



(91) 



For R(fi) < 2M < 2Mq eq. (l9TT) is real and thus the contribution (l90l) is real. If on the 
other hand f 2 % < fu we integrate p c i from fu to f 2 i + e keeping f 2 fixed. This time we 
recall that 

Pel = Pel + Pel (92) 

where p° cl is real because we are outside the interval (2M , 2H), and 



dT 

p cl = -{R'{h + e) - l)Z>(f 2 ) + (R(h) - f 2 )—V{f 2 ) 

or 1 



(93) 



with T now given by logt^. Again all the items in eq. f!93j) are real. Thus we can start 
e.g. with fu = T2% — £■ The integration along the path f\ = f 2 — e up to fo is very simple 
because 

1 + 



Pc2 +Pcl = f 



2M 



2H 



loe 



V2-P2 

f 



I — / 2M 



(94) 



whose complete discussion has already been done in Sect. El The "gap" is 2M, 2H. More 
difficult is the analysis for fo < f < 2H. Now we have 



with 



Pc2 + Pci = r 



2M 



P'l 



2H , 1 
— ; log — 

r 



H — Mq 



V2-P2 



(95) 



(96) 



Moreover 



AW (f 2 + e) + R'(f 2 + e)y/A* 

?2 



(97) 
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where 



and 



„ ML- M 
A = — 

r 2 



m 2 

2f| 



W(f 2 + e) 




1 



2M n 



(98) 



(99) 



r 2 r 2 

p[ diverges with positive residue at fx = 2H and thus also W(f 2 + e) and p' 2 of eq. flDTl) 
diverge like p[. As a consequence the analytic continuation of V 2 — p 2 below f = 2H is 
negative which makes the numerator of the argument of the logarithm in eq. ( 1951) negative. 
We show now that such numerator stays negative all the way for f < 2H. The argument 
of the square root in (1971) never vanishes so that at f 2 = 2M' Q the numerator in (1971) 
reduces to 

A(W(r 2 + e) + R'{r 2 + s)). (100) 

We can now explicitly compute W(f 2 + e) and R'(r 2 + e) at f 2 = 2M' Q . At such a point 
we have 



R\f 2 + e) = 1 



1 — / — 




(101) 



while 



W(f 2 +e) 



2H 1 H 



2M' Q 



r 2 I _ W. 



H 

K 



H 

A/,', 



1 




W I (102) 



which means that there is no pole in p' 2 at f = 2Mq. Thus j3 2 is regular below 2H and 
^2 — P2 cann ot change sign. 

In this way we proved that for the crossing of a null shell and a massive shell, into a null 
shell and another massive shell even with a change of mass, the imaginary part of the 
integral (181)]) is still given by eq. (18"9"]) . 

The reasoning in the general case of both mi and m 2 different from zero is dealt with 
similarly. It is sufficient to examine the case f 1 < f 2 the other case being now equivalent. 
We shall first discuss the integration along the path fi+e = f 2 up to tq. For ri+e = f 2 = f 
we have 



Pc2 + Pel 



2M 



2H 
— - log 

r 



\ — / AM. 



with 



H - M = V 2 + 



77*2 

2f 9 



P2 



2H 

f 2 



(103) 



(104) 



20 



which as before makes ]3 2 diverge at f 2 = 2H. For pi we have 



Pi AtWih + 0) + R'{h + G)^A{ - (1 

f, x — 2M ° 

fi 



2Mq \ iri\ 



(105) 



with 

M - M m| 



A = (106) 

p 2 diverges with positive residue at f = 2H and pi also diverges like p 2 as 

R\h + 0) = 1 + 4- (107) 
ri 

and W(f\ + 0) as given by eq. fllll) also diverges like p 2 . Then for f just below 2H we have 
both V 2 — ]3 2 < and Vi — pi < 0. As before all the point is in proving that pi is regular 
below 2H i.e. does not diverge at f = 2Mq. To this end we must examine the numerator 
of eq.CEUSD at f = 2M . We have for f = 2M 



W(2M + 0) = I( J^- - 1) + ^= > (108) 

2VM SM^^+l) 

as H > M . With regard to R'(fi + 0) which is negative for fi = 2if — £ it cannot change 
sign for 2Mq < f < 2if because at the point where R'{f\ + 0) vanishes 



2M 

W(f 1 + 0) = x /R' 2 (r 1 + Q)-l + -^ (109) 



would become imaginary while from W{f\ +0) = y^2H/fi+p2/ri we have that W is real. 
Then at f = 2M the numerator in eq. (11051) vanishes and pi is regular all the way below 
2H. Thus we are in the same situation as in the previously discussed case. Below f = 2M 
the argument of the logarithm in eq. (H03p becomes positive again due to the change in 
sign of the denominator 1 — \j2Mjr. The integration for f > fo is treated simply by 
exchanging m\ with m 2 . 
Finally we deal with the integral 

ffu+e 

Pc2 df 2 (110) 



r-2i 



which makes the two initial points coalesce. As in the previously discussed case of one 
massless shell, all the point is in proving that T> is real. In addition to the contribution 
of B which we already proved to be real, we have now the contribution of C. We have 
already proven that at r = f 2 — e 



R'(r) - jR'(ry -l + ^fl = R'( r ) - W(r) (111) 

R[r) 
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equals 




(112) 



which is negative for f 2 < 2H and thus in particular for f 2 < 2M. In the interval 
ri < r < f 2 < 2M the term fill II) cannot change sign because —1 + 2M /R is positive. 
Moreover we have 



R\h - e) - W(fi - e) = R\h + e) - W{f l + e) + — ^- - . (113) 



But we proved after eq. (11071) that below 2H, the analytic continuation of V\ — pi is 

negative, implying that (11131) is negative, like R'{r\ + e) — W{f\ + e). The outcome is 

that the discontinuity of C at f\, being given by the logarithm of a positive number, is 
real. This concludes the proof in the case of the emission of two massive shells. 

10 Conclusions 

The main issue of the present paper is the treatment of two intersecting shells of matter 
or radiation in general relativity, in a formalism, apt to compute the tunneling amplitude 
for the emission of two shells. In order to do so it is necessary to adopt a gauge which 
is more general than the one used P, [9] in the treatment of a single shell. In the usual 
treatments of the tunneling amplitude for a single shell, a limit gauge is adopted. Already 
at the level of a single shell we show that the complete action contains a term in which the 
mass of the remnant black hole plays a dynamical role. Such a term is unimportant if the 
variation of the action is taken with respect to the total mass of the system, keeping the 
mass of the remnant as a datum of the problem, but becomes essential if one varies the 
mass of the remnant keeping the total mass of the system fixed as done e.g. in [Sj. The 
reduced canonical momentum even in the single shell instance is gauge dependent but 
the tunneling probability turns out to be independent of such a choice. All the treatment 
is performed in the general massive case, the massless one being a special case. The 
tunneling results are independent of the mass of the shell. The adoption of a general 
non-limit gauge allows the extension of the formalism to two or more shells. In this 
instance both the intermediate mass and the total mass become dynamical variables in 
the sense that the reduced action contains terms proportional to the time derivative of 
them. We show how to derive the equations of motion of both the interior and exterior 
shell by varying the reduced action and this is done both in the massless and in the more 
complicated massive case. With regard to the computation of the tunneling amplitude 
it is possible to prove an integrability theorem which allows to deform the trajectory in 
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coordinate space to a contour which drastically simplifies the computation. Firstly one 
proves in such a way that the result is independent of the deformation defining the gauge 
introduced in Sect. [21 and secondly one finds in the general massive or massless case that 
the tunneling probability is given simply by the product of the tunneling probabilities for 
the independent emission of the two shells. Such a circumstance is interpreted p3] as the 
fact that in this model we have no information encoded in the radiation emitted by the 
black hole. 



Acknowledgments 

One of us (P.M.) is grateful to M. Porrati and V. S. Rychkov for stimulating discussions. 



Appendix A 

For completeness we report here some formulas which are useful in the text. The con- 
straints are given by [3J El E] 



RR" R' 2 
H t = —r- + — 4 



TC r = n R R' - 7i' L L -p 8(r-r), 



Ln 2 RR'L' n L ir R L 



— + \/p 2 L~ 2 + m? 8(r — f). 



L 2L 2R 2 L 2 R 
From the constraints it follows [3] that the quantity 

lA n 2 R R(R') 2 

M = — H — — 

2R 2 2L 2 

is constant in the regions of r where there are no sources as there 

M'- 



T * " RL Hr - 



The equations of motion for the gravitational field are [9] 



L = N 



R 2 



-f ) + (N r L)>, 



R 



Ntt, 



+ N r R', 



N 



R 2 



R' 



+ 1 + 



R 

2 p 2 8(r — f) 
L 3 \/p 2 L~ 2 + m 2 



N'R'R 



(114) 
(115) 

(116) 

(117) 

(118) 
(119) 
(120) 
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ir R = N 



Lti 2 l it L ir r 



R 3 



R 2 



R!_ 
L 



N'R 
~J7 



+ (N r 7r R y. 



(121) 



The equations of motion for f and p follow from the above equations for R, L, ir R and 7Tl 
and the constraints. In fact using the relation 



dR(f) 
dt 



R(r + e) + R'(f + e)? = R(r - e) + R'(f - e)t 



we have 



?A& + AR = 0. 
Using now the equation of motion (11191) and the constraints 



AR'(r) 



where V = \/p 2 + m 2 L 2 , one obtains 



V 
R ] 



N(f)p 



Air L (r) 



P 
L 



L(r)y/p 2 + m 2 L 2 (i 



N r (r). 



Using the relation 
dL(r) 



dt 



and 



L'(r + s)f + L(f + e) = L'(f - e)f + Lir — e) — Ur + L 
AR' AN' R(f) 



tAtxr = -N(r) 



L(r) L(f) 



N r (f)An R 



derived from (11211) and 



one obtains 



V 



p 



d L (r) dAn L 
-Att l - Lir) 



dt 



dt 



N(f)p 2 L' 



N> 



L 2 {r)^p 2 + m 2 L 2 {r) L(f) 



^p 2 + m 2 L 2 (f) +p (N r )'. 



122) 



(123) 



(124) 



(125) 



(126) 
(127) 
(128) 

(129) 



The occurrence of the average values U = [L'(f + e) + L'if — e)]/2 etc. in the previous 
equation, was pointed out and discussed at the level of the variational principle in [H] . In 
the massless case m = however, using again relations (I122fl26|127p one proves that the 
r.h.s. of eq. (II 291) has no discontinuity i.e. there is no need to take the average value in 
the r.h.s. of eq. (j!29p . In fact let us consider the discontinuity of the zero mass version of 
the r.h.s. of eq. ([1221) 

, V NL' V N> , ATrf . 



L 2 



+ N r 



(130) 
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being r\ the sign of p, i.e. 



p(^-^ + A^). (1 31) 



From eq. (11 181) and the equation of motion (j!25p we have 

AT AT ML 

V -AL> = -Air R - ^Att l - LAN r/ (132) 

and from eq. (11211) 

t]NAtt r = RAN' + NAR'. (133) 

Substituting the two into eq. (11311) and taking into account that Ahl = —p/L and AR' = 
—rjp/R we have that expression (I13ip vanishes. This fact was discussed at the level of 
the variational principle in [TO] . 

From the equations of motion (I118II119|120[ 11211) it follows that in the region where there 
are no sources 

i.e. combining with fl 1 1 7j) . we have that in the region free of sources M. is constant both 
in r and in t. 



Appendix B 

Equation of motion for f. 
1) Outer gauge 
In this case 



Pl R(AC-AB) r\-ar-:)-J^ + J™+log[l-J™)). (135) 



Using 



and 



we have 



OH V V2mJ \V V r 



, + *\i- x m ] (13 6) 



m> {f - e) = -wthr) = - v^ + ?l (137) 



' H- \ U~^ =ll-r4\ (138) 



OH V 2H V 2H V £ ^ fzH \V V f 

v V f 
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from which 



2) Inner gauge 
This time we have 



f = N(r m ) 



V V r 



p 2H 



and 



Using 



and 



we have 



i.e. 



H-M = V--- P J— 



, , /2# /2M , / 2H 

Pl = r \£(r + e)-\l— + d— -log 1- J — 



dM 



P 



p 2M 



-i 



dC 
~dR' 



ir + e) 



V2MfJ \V 
1 / I2M p 



W(f + e) 



dp l c 



dM V 2M 



} __ P_ / P_ _ 2M 
V \ V V f 



p I2M 



:i39i 



(140) 



141^ 



(142) 



(143) 



(144) 



(145) 



Appendix C 

In this Appendix we derive the equations of motion for the exterior and interior shell from 
the reduced action (|6"2"]1 for masses mi and m-i different from zero. 

First we consider the variation SH ^ and 5Mq = 0. Under such a variation the coefficient 
of the fi term is given, using 



9R -m f) dT - dR '-(R> i) dT 

dH- {R - ri) dH> dH- {R ' l) dH 



(146) 



by (W(f\ + e)W{f\ - e)) — multiplied by 

dH 



R' 
(R-h) 



W(h - e){R' - 1) + W{h +e)(R'-l + (R- h)^ + (R' - 



-W{h - e)R" + W{h + e){R" + ^R> + ^R") 



(147) 
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where R, R', R" stay for R(fi), R'ifi + e), R"(r 1 + e). From 



we find 



W{h - e) = W(ri + e) + 



Pi 



R(h) 
Pi 



'1481 



(149) 



dfl'(ri + e) R(h) 

which substituted into eg. (11471) makes it vanish. This is an expected result as the motion 
of the exterior shell must not depend on the dynamics which develops at smaller radiuses, 
but only on the two masses H and Mq. 
With regard to the terms proportional to f 2 in addition to 



a dp® 



r-2 



<:2 



dH 



we have the term given by 



multiplied by 



dT 



W(f 1 +e)W(f 1 -e) dH 



-W(r 1 -e)(R'-l) + 



(-(R'-l) + (R-h)§- 2 ) 
W(h + e)(R'-l+ ?j±(R - h) + - 1)) 



(150) 
(151) 

(152) 



(R - h) 



W(h - e)(^-(R' - 1) - R") + W{h + s)(^-(R' - 1) - R" + 
dr2 dr2 



dvi 
dR l 



i dT , . . dv\ , ,, dT , . 



where again R, R', R" stay for R(fi), R'ifi + e), R"{f\ + e). As a consequence of eq. (11491) 
the above expression vanishes. Adding the contribution of the boundary term we have 

d P ° c2 



r 2 



dH 



N{r r , 







(153) 



which is the single shell equation of motion ( 1261) . 

The equation of motion of the interior shell is obtained from the variation 5H = 0, 
8M 0. In this case we have no contribution form the boundary terms. For the 
coefficient of f\ we obtain 



1 



W(f l + e 
1 



+ R'R 



dT 



W{f t - e) \dM 



W{h + e) dM 
dvi 
~dR 



(R' - 1)+ 



(154) 



<9M n 



dM 



R" 



1 



W{h + e) W{fx-£) 



OR 



OR' 



27 



having used the relation (see eq. (|T3l) ) 



dF 



1 

if 



dMdR' DM 
Substituting in eq. (1154D the relation 

P 1 ^ + R W{H+e) 



(155) 



dvi 



dM R{R'{h + e)px - V 1 W(f 1 + e)) 
such a coefficient of f\ becomes simply 

1 



(156) 



f x R'{h + e) - W{f x + e) 



(157) 



which is the non zero mass generalization of the coefficient of f\ appearing in eq. (170l) . 
Then using eqs. ( |T5|29l) for N{f\) and N r (r{) we have 



= 

+r 2 



^R'(f + e) - W(f + e) 



-I -1 



dp 



d 2 F Vo 



dF 



dM OM dR'f 2 dM 



a: 



(158) 



^R'(f + s) 



+h dM- R{h + £ 



W(r + e) 
8 2 F 



Pi R(h + e) 
V x W{f + e) 



d 



dM dR' yn+£) T2 df 2 



dM, 



o 



where a = V 2Hf / (\/2Hr + p) . The first line is simply the equation of motion for f 1; the 
second line vanishes being simply the equation of motion for f 2 due to the variation of 
Mo (see eq. (l32|) ) while the sum of the third and fourth line vanishes in virtue of relation 

(EHD. 
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